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ABSTRACT 

We investigate the stability of incompressible, exact, non-ideal magnetorotational (MRI) modes against par- 
asitic instabilities for a broad range of dissipation coefficients, relevant to astrophysical and laboratory envi- 
ronments. This allows us to uncover the asymptotic behavior of the fastest parasites in terms of the Elsasser 
number. A,,, when viscous effects are not important. We calculate the fastest growing parasitic modes feeding 
off the MRI and show that they correspond to Kelvin-Helmholtz instabilities for A^, > 1 and tearing mode 
instabilities for A,, < 1. We study in detail the regime A,, ~ 1 where both types of modes present comparable 
growth rates. We conjecture about the asymptotic behavior of saturation based on the idea that the saturation 
level of the MRI can be estimated by comparing growth rates (or amplitudes) of primary and secondary modes. 
In the ideal magnetohydrodynamic (MHD) limit. A,, ^ 1, where Kelvin-Helmholtz modes dominate, these 
estimates lead to saturation levels for the stresses that are in rough agreement with current numerical simula- 
tions. For resistive MHD, A,, < 1, the stresses produced by the MRI primary modes, at the time when the 
fastest tearing modes have growth rates similar to their own, decay proportionally to the Elsasser number. This 
behavior seems consistent with numerical simulations of resistive MHD shearing boxes. 

Subject headings: accretion, accretion disks — black hole physics — instabilities — MHD — turbulence 



1. INTRODUCTION 

The transport of mass and angular momentum in ac- 
cretion disks remains one of the least understood pro- 
cesses i n modern astrophysics. The standard accretion disk 
mode l ("Shakura & Sunyaev" 119731: iLynden-Bell & Pringld 
ll974t [Frank et al. 2002) is based on the assumption that tur- 
bulence provides an efficient mechanism for enabling accre- 
tion but magnetic fields, thought to be crucial for driving the 
turbulence, do not play an explicit role. It is cur rently be- 
lieved that the magnetorotat i onal instability (MRI: | Velikhovl 
IT959I: iChandrasekha^ [19601: iBalbus & Hawlevl I1991L II 9*981) 
is responsible for driving the magnetohydrodynamic (MHD) 
turbulence required for efficient angular momentum transport 
in astrophysical disks. 

Since the appreciation of the relevance of the MRI 
to accretion disks, significant progress has been made 
in understanding the physics of the instability and in 
characterizing the ensuing turbulent state. A large set 
of numerical simulations (see, e.g., Hawlev et al] 19951: 
Brandenburg et al. 1995'; 'Miller & Stone 2000'; 'Fleming et al.' 
2000; Sano & Inutsuka 2001; Sano et al. 1998, '2004; 
Turner et al. 2007; Fromang et al. 20071: iLesur & LongarettH 
2007i: lObergaulinger et al., ,2009.) have provided insight into 



the turbulent MHD flows under various physical conditions. 
However, the fundamental processes that determine the 
strength of the turbulence in the non-linear regime are 
yet to be deciphered. At present there have only been a 
handful of studies addressing the theore tical aspe cts of this 
probl e m (see, e.g., [G oodman & Xu 199 4l:lKnobloc h & Julien 
20051: lUm urhan et al. 2007; Tatsuno & Dorland 2008; 
Jamroz et al . 2008a b; Latter etal. 2009; Vishniac 2009t 
Pessah & GoodmanI 12009) . Thus, understanding the mech- 



anisms that lead to the saturation of the MRI constitutes a 
fundamental problem in modern accretion physics. 

Building models that capture the essential physics of the 



MRI and its saturation is crucial in order to build angu- 
lar momentum transport models, and eventually global disk 
mod els, beyond the standard viscous accreti o n disk (see, 
e.g.. iKato & Yoshizawal [1991 lOgilvid [lool iPessah et alj 
I2006H) . The identification, an eventual understanding, of 
correlations and scaling laws born out of the synergy be- 
tween numerical simulations, analytical, and semi-analytical 
work ([P essah et al. 2006a, 2007, 2008; Blackman et al. 200^ 
iHubbard & Blackman 2008: Lesaffre et al. 2009.) is likely to 
provide important insight toward this goal. This could also 
prove to be a fruitful approach toward building sub-grid mod- 
els for related microphysical processes which are unfeasible 
to incorporate in g lobal simulations in a direct way (see, e.g., 
lArnett et alJl200S1) . 

In a recent paper, iPessah & GoodmarJ ( l2009l) provided a 
summary of a parametric study of incompressible, MRI- 
parasitic instabilities in a viscous, resistive plasma in dissi- 
pative regimes accessible to current numerical simulations. 
They stated that the fastest growing modes are related to 
Kelvin-Helmholtz and tearing mode instabilities. In this 
work, we provide the details of the solutions to the differ- 
ential equations involved and solve for the fastest growing 
modes in order to provide support for these assertions. We 
identify the existence of a critical Elsasser number of or- 
der unity and show that Kelvin-Helmholtz and tearing modes 
dominate in the quasi-ideal and resistive MHD regimes, re- 
spectively. By means of a systematic study of the param- 
eter space involved, we reveal scaling laws that describe 
the behavior of the fastest growing parasites in these lim- 
its. The joint analysis of the asymptotic behavior of the 
MRI and the parasitic modes that feed off it provides insight 
into the characteristics of a viable saturation mec hanism in 
regimes that ar e relevant to astrophysic a l (see, e.g. [jinll 19961: 

OarnmiJ 1 1 QQ/^l 'Qcinrv Xr \/r;,7cim J ' 1 QQC* lQci1moi-r,n Xr Wai-rllo 



Gammie"1996'; 'Sano & Mivama' '1 9991: ISalm eron & Wardle 
2005; Brandenburg & SubramaniaJ "[20051; [Balbus & Henrj 
20071) as well as experimental environments (see, e.g. IJi et al. 
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200ltlGoodman & Jil llOOl ISisan et al.ll2004t iLiu et al.ll200l 
Rudiger et alJl2003h . 



2. ASSUMPTIONS AND PRIMARY MRI MODES 

Consider a cylindrical background characterized by an an- 
gular velocity profile ft = il{r)z, threaded by a vertical 
magnetic field B = B^z. We work in the incompressible 
limit which is relevant when the magnetic fields involved 
are so weak that the saturation of the MRI occurs at sub- 
equipartition fields and consider the effects of constant kine- 
matic viscosity h' and resistivity rj. The equations governing 
the dynamics of this MHD fluid in the shearing box approxi- 
mation are given by 



dv 
'dt 



El 

Sir 



{B-V)B 

A'Kp 



i^V^B, 



(1) 
(2) 



— + KV)B = (B.V) 

where P is the pressure, p is the (constant) density, and the 
factor 

9 = - 



dlnr 



(3) 



parameterizes the magnitude of the local shear at the fiducial 
radius ro. The continuity equation reduces to V-v — and 
there is no need for an equation of state since the pressure can 
be determined from this condition. 

The set of Equations ([T]i and (|2]) possesses exact solutions 
of the form 

i; = -gr!o(r-ro)0+Fosin(ifz)e", (4) 
B = B,^z + BaCOs{Kz)e^\ (5) 
against perturbations parallel to the background magnetic 
field. For a given wavenumber K perpendicular to the disk 
mid-plane, the growth rate T satisfies the dispersion relation 

where T,, = T + j/A'^, r„ = 1 



^Tl)~AK^vl^nl = , (6) 

+ fjK'^, K = y/2{2~q)no 
z/y/Anp is the Alfven 



is the epicyclic frequency, vaz = 
speed, and Hq is the local Keplerian frequency. The ratio of 
the amplitudes Vq / Bo of the vectors characterizing the MRI 
velocity and magnetic fields, 

yo = Vb (cos6lv,sin6lv,0) , (7) 
Bo = So(cos0B,sin0B,O), (8) 
as well as their directions 6y and 9 b, see Figures [l] and [1 1 
are known functions of {v, t], K), see iPessah & Charil (l2008l) 
for a detailed discussion. The growth rate F has a unique 
maximum, T^^^{v, 77), at K Krm,y,{v, 77). 

Unless otherwise mentioned, we work with dimensionless 
variables defined in terms of the characteristic length and time 
scales set by the background Alfven speed and the local an- 
gular frequency: Lq = VAz/^o andTo = l/J7o- We subsume 
the scales related to viscosity and resistivity into the dimen- 
sionless quantities 



_ V\z 



^Az 



(9) 



(10) 



whose ratio is the magnetic Prandtl number, Pm = A^/ = 
ly/rj. The quantity A^^ is known as the Elsasser number, while 
its viscous counterpart A^, is related to the Reynolds number. 
In our dimensionless variables, magnetic field strengths are 
defined relative to the background field Bz, while A^, = v^^ 



and A^ = 



3. PARASITIC INSTABILITIES 



The exact equations for the dynamical evolution of the ve- 
locity and magnetic fields corresponding to secondary insta- 
bilities bv(x^ t) and SB{x, t) affecting an exact MRI mode 
are obtained by substituting in Equations ([T]i and (|2]i the 
ansatz 

V = -qno{r - ro)0 + Vq sm{Kz) e^* + Sv , (11) 
B = BzZ + BoCOs{Kz)e^* + 5B . (12) 

This procedure leads to a set of partial differential equations 
in space and time whose solution is beyond the scope of the 
present paper. In order to make the problem of the stabil- 
ity of the non-ideal MRI modes against parasitic instabili- 
ties more tractable we adopt the same set o f assumptions 
and approximations s tated in P essah & GoodmanI (I2009i) . see 
also iLatter et al.' ('2009), which are similar to those adopted 
by Goodman & Xu ( 1994) for ideal magnetohydrodynamics 
(MHD). In brief, we assume that the amplitudes of the pri- 
mary modes are large enough that we can neglect the influ- 
ence of the weak vertical background field, the Coriolis force, 
and the background shear flow on the dynamics of the sec- 
ondary modes. The source of free energy for the secondary 
instabilities should increase with the amplitudes of the pri- 
mary modes. Therefore, the secondary growth rates should 
eventually outrun the primary MRI modes. We thus assume 
that we can neglect the temporal variation of the sinusoidal 
velocity and magnetic fields of the MRI modes. 

3.1. Equations of Motion for the Parasites 

Under the assumptions stated above, the equations govern- 
ing the dynamics of the secondary instabilities are 

dtSv + (At;-V) 5v + {5v-V) Av = -V{5P + AB-SB) 
+ {AB-V)SB + {SB-V)AB + i'V^Sv, (13) 
dtSB + (AvV) SB + {6vV) AB = 



(AB-V) 6v + {6B-V) Av + ijV^SB , 



(14) 



where 'V-Sv = 'V-SB — 0, SP stands for the pressure pertur- 
bation, and we have defined the time-independent amplitude 
of an unstable MRI mode with wavenumber K as 



Av = Vo sm{Kz) , 
AB = Bo cosiKz) . 



(15) 
(16) 



We seek solutions to Equations ( fTsT i and ( fT4l ) of the form 

6v{x, t) = Svo{z) exp [st — ik-x] , (17) 
dB{x, t) = SBo{z) exp [st - ik-x] , (18) 

Here, fc = fch + kzZ, where fch = k^x + kyfj is a horizon- 
tal wavevector, see Figures [T] and |2] and kz is a parameter 
with < kz/K < 1/2. The eigenvalue s determines the na- 
ture of the temporal evolution of a given parasitic mode and 
it must be solved for together with the amplitudes dvo{z) and 
SBo{z), which are 27r/_ft'-periodic functions in z. 
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Fig. 1. — Three-dimensional representation of a primary MRI mode. 
Tlie perturbations over tlie background Keplerian velocity and magnetic 
fields are given by Av = Vo(cos ©V: siii^V: 0) sin(ii'2) and AB = 
i?o(cos Sb I sin Sbi 0) cos(ii'z), respectively. For a given primary MRI 
mode with wavenumber K, the orientation of the planes containing Av and 
AB. as well as the ratio of the amplitudes Vo/Bq, are functions of the vis- 
cosity and resistivity. The wavevector fci^ characterizes the horizontal wave- 
length of a given parasitic mode and 6 denotes the angle between this vector 
and the radial direction, see also Fig.|2] 



The set of six first order differential Equations (fTST i and ( fT4l i 
can in principle be solved for the secondary velocity and mag- 
netic fields (SP can be eliminated using V-^t; — 0) by requir- 
ing that the "boundary conditions" 

Sv{x + 2tt/Kz, t) = 5v{x, t) exp {2nik^/K) , (19) 
SB{x + 27t/Kz, t) = 5B{x, t) exp {2Trikz/K) , (20) 

be satisfied for all z. We follow an alternative procedure that 
leads to a set of higher order differential equations for Sv^ and 

SB,. 

Substituting expressions (fTTT l and ( fTSl l into Equations ( fTsT i 
and (fT4l l we obtain for the z-component of the momentum and 
induction equations 

sSv, — i{ki^-Av)Sv2 — —dzSP 

-i{k^-AB)5Bz+v{dl-kl)5vz, (21) 




Fig. 2. — Two-dimensional representation of the projection of a primary 
MRI mode onto the (f , plane. The vectors Vo = Vb(cos 0v> sin^v, 0) 
and Bo = Bo (cos 6b i sin Sb , 0) represent the projection of the MRI ve- 
locity and magnetic fields, with associated direction angles 9\j and 6^. The 
versors fcjj and fcp characterize, respectively, the directions parallel and per- 
pendicular to the horizontal wavevector fch of a given parasitic mode. 



We can eliminate the pressure perturbations 5P by taking the 
derivative of Equation ( |23] | and substituting for dzSP from 
Equation ( |2TI ). The result reads 

[s - i[k^-Av) - v{dl - kl)\{dl - k\:)Svz 
+ i{k^,-dlAv)6vz ~ i{k^,-dlAE)5Bz 



+ i(k^-AB)(dl-k\)bBz =0. 



(25) 



Using the explicit form of At; and AB from Equations ( fTST i 
and ( fT6l ), we obtain a complete set of equations for the z- 
components of the secondary amplitudes Sv, and SB, 



- iK'^{k^-Av)5vz + iK^{ki,-AB)SBz 
+ i{k^-AB){dl ~ kl)5Bz ^ , 



(26) 



[s - i{kh-Av) - ?7(af - kl)]6Bz + i{kh-AB)6vz = . 



(27) 



sSBz-i{kh-Av)5Bz 
-i{ki,-AB)Svz 



v{d!-kl)SBz 



(22) 



The horizontal component of the momentum equation reads 

sSvh~ i{kh-Av)6vh + {dzAv)Svz — —'VhSP 
-i{kh-AB)5Bh + {dzAB)5Bz 
+ vdlAv + v{dl ~ kl)ivx, , (23) 

where 6v\i = ki^-Sv and (5Bh ~ k^-SB. Taking the scalar 
product of this equation with ik^, and using the divergence- 
less nature of the secondary velocity and magnetic fields to 
write i{k\i-5v) — dzSvz and i{k\i-SB) = dzSBz, we obtain 



sdzSvz — i{kh-Av)dz6vz + i{kh-dzAv)5vz — 
- i{ki,-AB)dzSB, + i{k^-dzAB)SBz 
+ iydl{ki,-Av) + v{dl - kl)dz5vz ■ 



-kl5P 



(24) 



Finally, defining the differential operator 



(28) 



and regrouping terms, w e obtain the equations of m otion for 
the parasites presented in lPessah & GoodmanI (l2009t) : 



(s -I- vQ.)Q5vz - z(fch-At^)(Q - K'^)5vz 
+iikh-AB){Q-K^)5Bz = 0, 



(29) 



(s + riQ)SBz + i{ki,-AB)5vz - i{ki,-Av)5Bz = . (30) 

Thus, the set of assumptions stated in the beginning of this 
Section allows us to reduce the problem of analyzing the sta- 
bility of the exact MRI modes against secondary perturbations 
in terms of a set of linear, ordinary differential equations with 
periodic coefficients. 
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3.2. The Eigenvalue Problem and Its Solution 

The set of Equations ( |29] l and dSOl ) can in principle be inte- 
grated along the z-coordinate subject to the boundary condi- 
tions ( fT9] l and ( l20l i. Note, however, that the terms proportional 
to the viscosity v and resistivity -q constitute singular pertur- 
bations in the limit v,rj 0, i.e., they multiply the highest 
derivatives in the differential equations. It is thus convenient 
to work in Fourier space and seek solutions of the form 

oo 

5v,{z-k,)^ J2 a„e-^("^+'=^)^e-''^^-=^ , (31) 

n — — CX3 
oc 



Presuming that fch > 0, the operator Q is Hermitian and 
positive definite on functions f{z) that are sufficiently well 
behaved for z G (— oo, oo): specifically, / and its first deriva- 
tive are square-integrable. When / represents Svz or SBz, 
these mathematical conditions correspond to the physical re- 
quirement that the energy in the parasitic mode be finite; 



2 
1 



{\Sv\ 



\6Bndz. 



[\Svz\-' + \dBz\' + k-^i\djvz\' + \dJBz\^)] dz . 

(33) 

The series OTT l and (|32] | have infinite energy because they do 
not vanish as \z\ oo; however, superpositions of such so- 
lutions with a small range of can be made to vanish at 
infinity. 

When the operator Q acts on the Fourier Series each indi- 
vidual terms incurs a factor 



)„ = kl + {kz + nKf 



(34) 



Thus, the differential equations ( [29l ) and (l30b lead to recursion 
relations for the Fourier coefficients a„ and /3„ of the form 

= 2(s + vQn)Qnan 

- i (fch-Bo)[(Q„-l - K^)l3n-l + [Qn+l - K^)Pn+l] 

- (feh-Vo)[(g„-i - K^)an-i - {Qn+l - K^)an+i] , 

(35) 

= 2(s + r/Q„)/3„ 

- i {kh-Bo){an-i + - {kh-Vo)iPn-i - Pn+i) ■ 

(36) 

It is convenient to use the natural scales provided by the pri- 
mary MRI mode and rescale the variables; Qn Qn/K^, 
kh K/K, Un an/Bo, Pn Pu/Bq. We thus obtain 



KBn 



KBo'^" 
^ (fch--Bo) 

iK-Vo) Vq 

2Qn Bq 



[(Q„_l - l)/3„„i + [Qn+l - l)Pn 



+1. 



[{Qn-1 - l)ari-l - {Qn+l - l)a„+l 



Pn = --pTTrQnPn + « ^ (ctn-l + Oln+l) 



KBo 



KBo^ ' 2 

(fch-Vo) Vq 

2 Bo 



(/3„_i -/3„+i).(38) 



where the versors Vq = Vq/Vq and Bq = Bq/Bq provide 
the direction of the MRI velocity and magnetic fields ' . 

The system of coupled linear Equations ( [37| i and ( [38] ), with 
n = — oo, . . . , oo, can be written in matrix form as Mq ~ sq, 
where M is a band-diagonal, complex, non-Hermitian ma- 
trix and the components of the eigenvector q are defined ac- 
cording to Qn = an/2 for u ^ 2m and g„ = /3„/2+i for 
n = 2m + 1. Boundary conditions such that (q;„,/3„) 
must be imposed on the eigenvectors in order to ensure that 
the Fourier Series ( 1311 1 and ( |32] | converge as \n\ oo. 

In practice we set q;n = /3n = for |iV| > A^max and di- 
agonalize the 2(27Vi„ax + 1) x 2(27Vmax + 1) matrix M for 
a given set of values (i/, rj, K, Bo,kz,kii,9). The core of the 
algorithm employed to diagonalize the matrix AI is based on 
the LAPACK routine ZGEEVX, using the option of apply- 
ing a balancing transformation to improve the conditioning 
of the eigenvalues and eigenvectors. We solve for the finite 
set of 2{2Njnax + 1) eigenvalues {s} and eigenvectors {q} 
(if needed), taking increasingly large values of A^max until 
convergence to the desired accuracy is reached; Ny^ax = 30 
seems to do a very good job across the parameter space that 
we explored. 

3.3. Parameter Space Survey 

We consider a grid defined by the parameters = 
{1, 10, ... , 10'^} andPm = {10"'^, 10"^ . . . , 10^} and solve 
the system of Equations ( l37l i and ( |38] | as a function of the 
variables (Sq, fc^, 0, fch) focusing our attention on the stabil- 
ity of the fastest growing primary MRI mode ifinaxl*^, v)- For 
each pair of values (A^, Pm), we solve for s(i?o, ^z, 6, kh) 
over a grid of values according to kz = {0.0, 0.1, ... , 0.5}, 
9 = {0°, 5°, . . . , 180°}, and = {0.0, 0.01, . . . , 1.0}. We 
find the fastest growing mode for a fixed value of the primary 
MRI magnetic field and iterate over Bq, in steps of 0.1, until 
the fastest secondary instabilities matches the growth rate of 
the fastest MRI mode, i.e., s,nax(t^, ?7, ^max) = TmaxCj^, ??)■ 
We denominate the corresponding value of Bq by i3g^'(i/, 77) 
and we refer t o it as the " saturation" amplitude of the MRI 
magnetic field dPessah & Go odman 2009)^. 



The real and imaginary parts corresponding to the eigenval- 
ues with fastest growth rates, Sjnax{'^, ?/, -f^^max, kz, 0), of var- 
ious secondary modes that feed off the fastest primary MRI 
mode are shown in the upper and lower panels of Figure [3] 
respectively. The different panels in this figure correspond to 
Pm = {10^4, 10-3, 10-2, 10-1} with Aj, = 10^, and thus 

' Recall that the ratio of the velocity to magnetic field amplitudes, Vo /Bq, 
for the primary MRI mode is not an independent variable; it can be calculated 
for each set of values {i/, rj, K), see Eq. 54 in Pessah & Chan (2008). It is 
convenient to scale the Fourier coefficients using Bq, rather than Vq, because 
the ratio Vo / Bq is proportional to A,, for A,, "C 1, see Eq. (49). 

^ A criterion defined in terms of matching amplitudes of primary and sec- 
ondary modes leads to MRI fields that are factors that are ~ 3 to ~ 10 times 
larger than _Bq^'(;/, r]). Note that in neglecting the temporal dependence of 
the background, it is assumed that the secondary growth rates are large com- 
pared to the primary growth rate; thus our definition of saturation implies an 
(37) extrapolation to the regime where this assumption is not strictly satisfied. 
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Fig. 3. — Real and imaginary parts for the parasitic instabilities with fastest growth rates normalized by the growth rate of the fastest MRI mode as a function 
of the oiientation of the horizontal wavevector fc^ with respect to the radial (8 = 0) direction. In all panels, the viscous counterpart of the Elsasser number is 
Av = 10^ ; the magnetic Prandtl number varies according to Pm = {lO"*, 10"^, 10"^, 10"^}, and thus the Elsasser number is A,, = {10"\ 1, 10, 10^}, 
from left to right. The various lines in each panel correspond to different values of the parameter kz. The MRI magnetic field Bq = Bq^^{i/, tj) is such that 
the fastest parasitic growth rate, maximized over feji, 6, and kz, matches the growth rate of the fastest primary !MR.I mode, rmax(i^; 7 7). The vertical solid lines 
show the angles 6\{u^r}) and 0-Q{y^r}) associated with the velocity and magnetic fields of the fastest MRI mode, see Figs. [T| and [2] The MRI velocity and 
magnetic fields are close to orthogonal, except in the case = 1; the vertical dashes fines show the directions perpendicular to 0b (left) and B^J (right). The 
parasites present fastest growth along the directions 6max ~ ^'v ™d ^max ~ Cb and they are associated with Kelvin-Helmholtz and tearing mode instabilities, 
respectively. The former type dominate for A,, > 1 while the latter prevail for A,, < 1. This is a general result that holds as long as A^/ > 1, see Section|4] 
Tearing modes are enabled by resistive reconnection and are thus expected to become increasingly relevant as A,, diminishes. The reduction of the growth rate 
of the Kelvin-Helmholtz instabilities at low Elsasser number is due to the decrease of the MRI velocity field ampfitude; Vb / Ba ~ A^ for A^ <g; 1. 



= {10"\ 1, 10\ lO^}, from left to right. The various 
curves in each panel correspond to different values of the pa- 
rameter kz; and Bq — Bf^^{v^ rj) in all of the cases shown. 

Figure [3] illustr a tes a generic result, described in 
iPessah & GoodmanI (12009 ') for a small range of {A^, Pm}. 
The fastest parasitic modes are non-axisymmetric (6'inax 7^ 
0), have the same vertical periodicity as the primary mode 
(kz = 0), and have purely real growth rates. The fastest 
modes have horizontal wavevectors, fch,inax, that are nearly 
aligned with either the velocity or the magnetic field of the 
primary mode. The first type dominate for A^ > 1 and corre- 
spond to Kelvin-Helmholtz modes that reach their maximum 
growth rates for 6'niax — ^'v. i-S-. the direction of the MRI 
velocity field, AV. The second type dominate for A^ < 1 
and are related to tearing modes that grow the fastest along 
the direction 6'inax — Ob, i.e., the direction of the MRI mag- 
netic field, A B. Note tha t because the angle 6'v — 6*8 — 7i'/2 
for A^ ^ 1 (iPessah & Chan 2008). the directions parallel to 
9\r and perpendicular to 9b do not differ by more than a few 
degrees. 

Table [T] provides information concerning the growth rates 
and geometric structure of both the fastest growing primary 
MRI mode and the associated fastest parasitic mode for a 



10^^. Tearing modes and Kelvin-Helmholtz modes 



range of Elsasser numbers A,, ~ {10 



10^}, with 



Ai, = 10"^. As shown in Figure |3] there is a sharp transi- 
tion in behavior around a critical Elsasser number of order 
unity. The low and high Elsasser number regime is charac- 
terized by asymptotic behaviors that are akeady evident for 



present the fastest growth at low and high Elsasser numbers, 
respectively. The physical properties of both types of modes 
are insensitive^* to the value of A^ provided A^, > 10. 

The manifest asymptotic behavior inferred from Table [T] 
simplifies the characterization of the physical nature of the 
parasitic modes as a function of the Elsasser number. Thus, 
we first analyze the asymptotic regimes al large and low El- 
sasser number in Section |4] and postpone the explicit analysis 
of the physical structure of the eigenmodes until Section|5] It 
will then be evident that it is only necessary to understand the 
nature of the modes with moderate Elsasser number A^ ~ 1 
and that the characteristics of the modes for A^ <C 1 and 
Ajj ^ I can be obtained using the scaUng relations presented 
in Section!!] 



4. ASYMPTOTIC BEHAVIORS AND SCALING LAWS 

Identifying in a systematic way the fastest growing parasitic 
instabilities as a function of the dissipation coefficients allows 
us to investigate their asymptotic behavior in various regimes. 
This in turn enable us to uncover simple scaling relations for 
their growth rates in terms of the amplitude of the MRI ve- 
locity and magnetic fields. In order to highlight the physics 
underlying the various scalings we work with variables with 

' The regime where viscosity plays an important role, i.e., A^ < 10, was 
explored in Pessah & Goodman 1 2009). Here, we mainly focus our attention 
on the regime where A^ > 10'^, which exhibits asymptotic behaviors that are 
more likely to be relevant to astrophysical and experimental environments. 
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TABLE 1 

Properties of Fastest MRI and Parasitic Modes 







^max 


9v[°] ' 


Sen 


^max [ ] 


^h,max/.^max 




Vo/Bo 


Type 


l.E-03 


7.50E-04 


8.66E-04 


0.0 


90.1 


90.0 


0.48 


4.9 


0.0017 


TM 


l.E-02 


7.50E-03 


8.66E-03 


0.0 


90.6 


90.0 


0.48 


4.9 


0.0173 


TM 


l.E-01 


7.39E-02 


8.55E-02 


3.5 


95.6 


95.0 


0.48 


4.9 


0.1689 


TM 


l.E+00 


4.28E-01 


5.16E-01 


22.1 


119.7 


30.0 


0.59 


4.5 


0.7187 


KH 


l.E+01 


6.96E-01 


8.88E-01 


40.6 


132.9 


40.0 


0.59 


3.8 


0.7816 


KH 


l.E+02 


7.44E-01 


9.59E-01 


44.5 


134.7 


45.0 


0.59 


3.8 


0.7751 


KH 


l.E+03 


7.49E-01 


9.67E-01 


45.0 


135.0 


45.0 


0.59 


3.8 


0.7743 


KH 


l.E+04 


7.50E-01 


9.67E-01 


45.0 


135.0 


45.0 


0.59 


3.8 


0.7742 


KH 


Note. 


— Data corresponding 


to the fastest g 


jrowing MRI mode and their associated fastest 


parasitic modes for Ai, 


= 10^. 



Note the sharp transition around A^ = 1 and the asymptotic behaviors at low and high Elsasser numbers; corresponding to 
resistive and ideal MHD and dominated by tearing modes (TM) and Kelvin-Helmholtz (KH), respectively. The asymptotic 
behavior for the modes is already evident for A,, ~ 10^^. 



physical dimensions in this section. 

4. 1 . Asymptotic Behavior of Primary MRI Modes 

In order to better appreciate the asymptotic behavior of the 
parasitic instabilities it is necessary to consider the asymptotic 
behavior of the primary MRI modes. The quantities that de- 
fine the structure of these modes, as well as some relationships 
between various timescales and amplitudes, are of particular 
importance and we briefly summarize them here. 



4.1.1. Ideal MHD 

In the case of a Keplerian shear profile the maximum 
growth achieved by the MRI in the ideal MHD regime, i.e., 

A^, > 1, is 



-no. 

4 ° ' 



(39) 



(40) 



The fastest growth rate is related to the Alfven frequency, 
uJAz = K^jiaxVAz, associated with the background magnetic 
field Bz, via 



and corresponds to the mode with wavenumber 
K 




with LUAz 



(41) 



corresponds to the mode with wavenumber 
K -,f^-A^ 

V 4 ' VAz 



(45) 



In this case, the various relevant inverse timescales are related 

via 



^UJAz 



where 



LOAz = ^maxI'Az = 

is the Alfven frequency and 




(46) 



(47) 



(48) 



is the inverse of the resistive timescale across a lengthscale of 
the order of the MRI wavelength, i.e., if^ax- The magnetic 
field and velocity field amplitudes for the fastest primary MRI 
mode are related via (see Eq. 92 in P essah & Chan.2008 ) 



Bn 



:A. 



(49) 



and thus Vq/Bq ^ in the limit A,, ^ 1. The planes con- 
taining the fastest growing MRI velocity and magnetic fields 
are characterized by the angles, see Figure |6] 



The magnetic field and velocity fiel d amplitudes for the fastest 
primary MRI mode are related via (iPessah et al.ll2006"ah 



(42) 




The angles characterizing the planes which contain the veloc- 
ity and magnetic fields of the fastest growing MRI mode are 
given by, see Figure |6l 



sm 



' cos Ob 



1 

71' 



(43) 



4.1.2. Inviscid, Resistive MHD 



In the inviscid, resistive MHD limit, i.e., A^ 1 and A^ <C 
1, the fastest MRI growth rate 



sin6lv = -A^, (50) 
cos6'b = -A^. (51) 

4.2. Asymptotic Behavior of Parasitic Growth Rates 



The dependence of the (dimensionless) growth rate of the 
fastest parasitic modes on the amplitude of the primary MRI 
magnetic field, Bq, is shown in Figure |4] The various sets 
of data points corresponds to all the pairs of values (Ai,. Pm) 
that can be formed with = {10^, • ■ ■ , 10''} and Pm = 
{10^^, 10^^, . . . , 10^}. The various colors (and sizes) corre- 
spond to different values of A,y. The values of A^ associated 
with each set of data points are shown on the right. Open cir- 
cles identify Kelvin-Helmholtz modes that feed off the MRI 
velocity field at 6* ~ 9y, which exhibit the fastest growing 
rates for A,, > 1. The growth rates for the Kelvin-Helmholtz 
modes is almost insensitive to dissipation for A^,A^ > 1. 

(44) 

Filled circles correspond to tearing modes, that feed off the 
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MRI currents at ~ 9b- These are the fastest growing sec- 
ondary instabilities for < 1. 

The following (dimensionless) equations, shown as dashed 
lines in Figure |4] provide an excellent description of the 
fastest growth rates associated with Kelvin-Helmholtz modes 

s^ax = O.2OB0 for 9~ev,A^>l, (52) 

Smax^O.lOBo for 9^ev,A^^l, (53) 

while the scaling of the growth rates of the tearing-modes is 

well described by 



O.ISBqA,, for 



A„ < 1 



(54) 



Although these are not shown in Figure |4] the growth rates 
of the (sub-dominant) Kelvin-Helmholtz modes obeys 



; = 0.37SoA2 for 9 



A„ < 1 , 



(55) 



There does not seem to be simple scaling relations for the 
(sub-dominant) tearing modes in the regime ^ 1 and 
A^ > 1 fore* ~ ^B. 

The critical value of the Elsasser number at which the 
fastest growing tearing modes (along the direction 9 ~ ^b) 
grow faster than the fastest growing Kelvin-Helmholtz modes 
(along the direction 9 ~ 6'v) is very close to (but slightly less 
than) unity, see Figure [3] i.e.. 



A' < 1 . 



(56) 



The existence of a an order unity critical Elsasser number that 
dictates the linear evolution o f the MRI has already been ap- 
preciated on both numerical ( Sano et al.lll99 8't Fleming et alj 
2000) and analytical grounds (see, e.g., ISano & Mivama 
1999t iPessah & Chad ifoOSh . Here we posit that this criti- 



cal Elsasser number also distinguishes which type of parasitic 
instability dominates the subsequent evolution of the MRI. 
In the remainder of this Section we argue, and demonstrate 
in Section |5] that there are two regimes such that the fastest 
growing parasitic modes correspond to 

A,, > A^ ~ 1 Kelvin-Helmholtz , (57) 

A,, < A^ ~ 1 ^ Tearing Modes . (58) 

4.2. 1 . Ideal Kelvin-Helmholtz Modes 

For the fastest growing Kelvin-Helmholtz instabilities with 
A^ > 1 and A^ ^ 1, the ratio of the horizontal parasitic 
wavenumber to the wavenumber corresponding to the fastest 
growing primary MRI mode depends very weakly on either 
Ajy or A,,, see Table [T] 

~ 0.59 for A, » 1, A^ > 1 . (59) 

max 

Therefore, restoring the physical dimensions into Equa- 
tion ( |52] | and using Equation (l42l) leads to 

s,„ax ^ 0.44 '^^EMlh. ^ 0.45 c^KH.o , (60) 

where ujaz is given by Equation dTTT l and 

t^KH.O = fch Vo . (61) 

Note that tiJKH.o corresponds to the growth rate associated 
with a Kelvin-Helmholtz instability feeding off a velocity 
field discontinuity of amplitude Vq. The amplitude of the MRI 
magnetic field needed in order for the growth rate of the para- 
sites to match the growth of the primary mode can be obtained 
by equating ( [39l ) and ( l52b . which leads to 

Bi5^*~3.8 for A^>1,A^>1, (62) 

in agreement with the results of lPessah & Goodmanl(l2009l) . 



10^ 



10" 



10" 



10" 



s -# ® • ' 



O Kelvin-Helmholtz 
# Teoring Modes 
- - Anolyticol Expressions 



A, >1.00 
A, =1.00 

A, =1.0e-01 
A, =1.0e-02 
A, =1.0e-03 
A, =1.0e-04 
A, =1.0e-05 



• A„=1.0e+07 

• A„=1.0e+06 

• A^=1.0e+05 

• A„=1.0e+04 
A„=1.0e+03 

• A„=1.0e+02 
1 1 1 1 L_ 



10 



B„ 



Fig. 4. — Growth rates for the fastest growing parasitic instabilities as a 
function of the amplitude of the MRI magnetic field Bg. The data shown 
corresponds to all the pairs of values (A^,Pm) that can be formed with 
= {102, . . . , lO''} and Pm = {10"'', 10-"^, . . . , lO''}. The various 
colors correspond to different values of A,/ for fixed Pm. For the sake of 
clarity, the size of the data points also decreases with Ai, for fixed Pm. The 
values of A,, associated with each set of data points is shown on the right. 
Open circles represent Kelvin-Helmholtz modes, which exhibit the fastest 
growing rates for A,j > 1. The growth rates for the Kelvin-Helmholtz modes 
is almost insensitive to dissipation for A,^, A,, > 1. Filled circles represent 
tearing modes, which exhibit the fastest growing rates for A,, < 1. In all 
cases, the data points are drawn for Bq < Bq^' (i^, 77). The dashed lines, 
proportional to Bg , represent fits to the data whose analytical expressions are 
provided in Section|4] 

4.2.2. Non-ideal Kelvin-Helmholtz Modes 

In the regime A^, 1 and A,, < 1, the fastest growing 
Kelvin-Helmholtz instabilities grow slower than the fastest 
growing tearing modes at the same MRI field amplitude. The 
ratio of the horizontal parasitic wavenumber to the wavenum- 
ber corresponding to the fastest growing primary MRI mode 
depends very weakly on either A^, or A,,, see Table[Tl 



0.59 for At, > 1, A,, < 1 . 



(63) 



Therefore, restoring the physical dimensions into Equa- 
tion (I53I 1 and using Equation ( l49l l leads to 



^0.37 



^KH,0 WAz 



0.42tJKH,0- 



(64) 



where ujkh.o is the Kelvin-Helmholtz growth rate defined in 
Equation i6l\ and ljaz is the Alfven frequency defined in 
Equation (|47] |. Note that despite the fact that these are Kelvin- 
Helmholtz modes, the resistive timescale enters because the 
resistivity implicitly affects the amplitude of the MRI veloc- 
ity field which is the energy source for these parasitic modes. 

4.2.3. Tearing Modes 

For the fastest growing tearing modes with A^, ^ 1 and 
A^ < 1, the ratio of the horizontal parasitic wavenumber to 
the wavenumber corresponding to the fastest growing primary 
MRI mode depends very weakly on either A^, or A^, see Ta- 
ble[Tl 



0.48 for A„>1,A^<1. 



(65) 
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A^=1.0e+07 
A,= 1.0e+06 
A.= 1.0e+05 
A„=1.0e+04 



A,= 1.0e+03 
A„=1.0e+02 
A„=1.0e+01 
A„=1.0e+00 



10" 



10" 
Pm 



10= 



Fig. 5. — Magnetic energy density corresponding to the fastest growing pri- 
mary MRI mode wlien tlie growtli rate of tlie fastest parasitic mode matclies 
rmaxCf, i.e.. Bo = Bq^^ . The data shown coiTesponds to all the pairs 
of values (Ai/, Pm) that can be formed with A,^ = {1, 10, . . . , 10^} and 
Pm = {10~^, 10~^, . . . , 10^}. The various colors correspond to differ- 
ent values of A^. There are two clear asymptotic regimes that con'espond 
to A,, > 1 and A,, < 1, which are associated with Kelvin-Helmholtz and 
tearing modes respectively, see Sectiong] For A^ = {1, 10}, the viscous 
quenching of the Kelvin-Helmholtz modes leads to larger saturation ampli- 
tudes for the primary MRI modes; note that these numbers correspond to a 
highly viscous regime. 

Therefore, restoring the physical dimensions into Equa- 
tion ( |54] | leads to 



iO.Sl 



LU,-, 



0.36 WAo , 



where we have used Equation (l46T l and defined 

khBo 

^AO ~ n ■ 



(66) 



(67) 



as the Alfven frequency associated with the horizontal MRI 
magnetic field Bq. 

The amplitude of the magnetic field that the MRI needs to 
grow in order for the growth rate of the parasites to be as large 
as that corresponding to the primary mode can be obtained by 
equating (l44l i and ( l54l i. We obtain 

5.0 for Ai.>l,A^<l (68) 



Dsat 
^0 



We can formulate a heuristic argument for the existence 
of an asymptotic limit for the amplitude Bq in this limit. 
The analysis of the classical tearing mode instability, where 
Bq ~ Bota.Bh{Kz), with K fixed, leads to growth rates 
7 ~ {'qK^)"{KBa)'^-" with a = 1/2 or a = 3/5 depend- 
ing on the details of how the problem is solved in the resis- 
tive layer For r] ^ 1 ^ v, there are only two characteristic 
scales in the problem, the resistive timescale [rjK^)^^ and the 
Alfven timescale (if i?o)~^- Therefore, the growth rate of the 
pai-asites must be s ~ {t]K^)°'{KBo)^-°' = 7;"ifi+"B^"", 
with < a < 1. At a fixed scale K, the growth rate of the sec- 
ondaries is proportional to a positive power of the resistivity, 
s ^ 77". However, if the scale K is given by the fastest grow- 
ing MRI mode then K ~ rj^^ and thus Smax ^ (^0)^ "/?/■ 
In this case, because Fmax ^ V^^^ the amplitude Bq at which 
Smax = Tinax must be independent of rj, i.e.. 



4 ^ 



l-Q 



(69) 




Pm 

Fig. 6. — Plot of the function | sin 0b cos 6b\, where 6b coiTesponds to the 
angle subtended by the magnetic field component of the fastest primary MRI 
mode and the radial direction. The data shown con'esponds to all the pairs 
of values (Ai/,Pm) that can be formed with Ai, = {1, 10, . . . , 10^} and 
Pm = {10~^, 10~®, . . . , 10^}. The various symbols (colors) correspond 
to different values of Ai, . In the ideal hmit, Ai, , A,, 2> 1, | sin 9b cos 6b I ~* 
l/\/2, while in the inviscid, resistive limit, i.e., A^ 2> 1 and A^ <C 1, 
I sin 9b cos 6b \ — > A,, . 

Therefore, Bq const, for A^ ^ 1. 

Because the growth rate of the parasitic modes and the 
MRI decay at the same rate for A^ <C 1, the amplitude to 
which the primary MRI magnetic field needs to grow to in 
order for both growth rates to become comparable is inde- 
pendent of the Elsasser number. This means that, as resistiv- 
ity increases the magnetic field that can be generated by the 
MRI before the tearing modes become dynamically important 
reaches an asymptotic value. In this respect, one could argue 
that the ability of the MRI to tap energy from the differential 
rotation does not inherently diminish as resistivity increases. 
The instability can achieve this goal via unstable modes with 
increasingly longer wavelengths that operate at increasingly 
slower rates. The situation is quite different for the stresses, 
however. 

4.3. Asymptotic Behavior of MRI-Stresses 

We calculate the dimensionless stress at saturation as 
asat = T^,fl{L,nQ)\ where f-* = i?^^' - is the 

sum of the Reynolds and Maxwell stresses 

^"^''^ VC{z)V^:l{z)dz , (70) 



1 



psat 



Ml 



1 



Bl%z)Bl%z)dz 



(71) 



Integrating these expressions we obtain the dimensionless 
stress dsat in terms of the parameter"^ (3 = ^'^Ll/v\^, 

"sat ^^[{Vt'f sin^v cos^v - {Bl^'f sin^B cos 0b] • 

(72) 



* IPessah & Goodmari j2009l) defined Osat = T^f/{SL^)^ and 

13 = S^L^/jj-^,, with S = 3Q0/2, in order to compare results with 
ILesur & Longaretg <2007ft . The definitions adopted here do not affect the 
expression or the numerical value for the product a/3. 
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Fig. 7. — Predicted values for the product «sat/3sat if saturation occurs 
when the fastest parasitic and primary MRI growth rates match. The data 
shown corresponds to all the pairs of values (A^ , Pm) that can be formed 
with = {1,10, ...,10^} and Pm = {lO"'^, 10"'', . . . , lO'^}. The 
vaiious symbols (colors) coiTespond to different values of A,^. In the ideal 
limit, A^, A^ 2> 1, Osatftat — * 0.4, while in the inviscid, resistive limit, 
i.e., Ai, 2> 1 and A,, 2> 1, Cfsat/3sat — > 0.5A,,. Despite the fact that the 
magnetic field at saturation asymptotes to a constant value, see Fig. |5] due 
to the behavior of the function | sin 9b cos 9b I for ^j? ^ 1. see Fig. |6] the 
dimensionless stress decreases at low Elsasser number. 

Using Equations (|43]l, and for the ideal MHD limit 
and Equations ( |49] l. ( BTl i. and ( |68] | for the inviscid, resistive 
MHD limit, we obtain the asymptotic limit of the dimen- 
sionless stress at saturation in Equation |72l as asat/3 — 6 for 
A,, ^ 1 and asat/3 — 12 A,, for A,, <C 1. We can relate the 
initial value of the /^-parameter, (3 oc l/Bf, with an estimate 
for the corresponding value at saturation via /3sat = P/Bq. 
Using Equation (|62] | and (l68l l this leads to 

asat/3sat=i0.4 for A^ > 1 , (73) 
asat/3sat=i0.5A^ for A,,<1. (74) 

It is worth mentioning that numerical simulations of MRI 
driven turbulence (with no explicit dissipation) carried out 
over a wide range of physical conditions lead to saturation 
values for the parameters gga t and 0sat that vary by several 
orders of magnitude (see, e.g. iPessah et al.ll2006al and refer- 
ences therein). However, the ir product r emains rough ly con- 
stant with Qsat/3sat ^ 0-5 dHawlev et al. 1995; San o et all 
l2004tlBlackman et al.ll2008h . 

We emphasize that, in our calculations, it is the behavior 
of the angle 9b 7r/2, see Figure |6] and not of the satura- 
tion amplitude _Bq** — > const., see Figure |5] as the Elsasser 
number decreases that brings about the asymptotic decay of 
the stress Tr^, oc A,, for A,, ^ 1, see Figure |7] (Because 
Vo/Bq ~ A,,, the contribution of the Reynolds stress to the 
total stress is negligible in the very resistive limit.) 

5. KELVIN-HELMHOLTZ AND TEARING EIGENMODES 

Throughout this paper we have stated that the fastest par- 
asites are related to Kelvin-Helmholtz and tearing mode in- 
stabilities in the regimes A^ > 1 and A^ < 1, respectively. 
The reason to delay presenting the rigorous evidence support- 
ing these claims up to this point is based on the existence of 
the asymptotic regimes and scaling laws presented in Sec- 
tion |4] We can now focus on the region of parameter space 



with "moderate" Elsasser number, A^ ~ 1, where changes 
in A,, produce non-trivial modifications to the parasitic mode 
structure. With this knowledge, and the insight gained in Sec- 
tion m it is straightforward to describe the structure of these 
modes in the limits of large and small Elsasser numbers. 

As a prelude to the calculations and description of the 
eigenmodes associated with Equations ( |37] | and dJSl l. we 
briefly describe some of the key physical signatures that can 
be used to recognize Kelvin-Helmholtz (shear driven) and 
tearing mode (resistive) instabilities. 

5.1. Signatures of Kelvin-Helmholtz and Tearing Mode 
Instabilities 

The simplest example of a physical system subject to the 
Kelvin-Helmholtz instability is that of a uniform, incompress- 
ible, shear flow with a discontinuous velocity profile V — 
UqX for a; > and V = —UqX for x < 0. This flow can 
be envisaged as a vortex sheet with an associated vorticity 
field u) — u!oy for z = and uj = for z 0. The en- 
suing Kelvin-Helmholtz instability emerges as a perturbation 
in the vorticity field. For a given wavenumber, the instabil- 
ity can be visualized as an alternating set of perturbations, 5uj 
which tend to incre ase/decreas e the local background vortic- 
ity (see Fig. 7.1.2 in lBatcheloili 2000). The Kelvin-Helmholtz 
instability reaches the largest growth rates in the absence of 
dissipation. In the ideal limit, the instability grows at a rate 
proportional to the magnitude of the wavenumber parallel to 
the direction of the velocity field and the magnitude of the 
velocity discontinuity. Finite viscosity modifies this picture 
in a thin boundary layer but the unstable modes with longer 
wavelengths are unaffected. 

The tearing-mode instability belongs to the class of sponta- 
neous (as opposed to driven) resistive instabilities. These are 
instabilities for which non-zero resistivity allows the plasma 
to access configurations with lower energies by enabling mod- 
ifications to the field structure (i.e., magnetic reconnection) 
that would be topologically inaccessible due to ideal con- 
straints. Perhaps the simplest physical system subject to the 
tearing mode instability is that of a uniform, incompressible, 
plasma threaded by a magnetic field that reverses direction in 
a narrow region i.e., a current sheet, and is uniform far away. 
The magnetic field given by B cx ta.nh{z / l)x is usually em- 
ployed as a convenient model to understand the physics in- 
volved in the stability of a thin current sheet j = jo{z)y 
of characteristic width I. The detailed solution of this prob- 
lem depends on the careful analysis of the physical pro- 
cesse s that take place in the reconnection region ( Furt h et al.l 
Il963h . However, asymptotic analytical treatments that match 
approximate solutions inside and outside the narrow region 
where the magnetic field reverses provide importa nt insight 
into the generic features of the solutions (see, e.g.. ISturrockl 
Il994t iBoyd & SandersonI 12003 ). A robust result is that the 
timescale for the growth of the instability is intermediate be- 
tween the Alfven timescale and the (usually much slower) re- 
sistive timescale across the reconnection region. 

The tearing mode instability is driven by the Lorentz force 
that acts toward the resonant (or null) surface where the mag- 
netic field goes through zero, i.e., z — 0. Non-zero resistivity 
allows the field lines to slip through the plasma enabling the 
reconnection of field lines of opposite polarity. This process 
takes place in various points along the null surface (determin- 
ing the wavelength of the unstable mode) and gives rise to 
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the formation of magnetic islands. These are a set of nested 
magnetic surfaces whose magnetic axis lie along the direction 
parallel to the background current sheet and are centered in 
the so-called "0-points". The magnetic islands are confined 
by a separatrix whose intersections with the null surface de- 
fines the so-called "X-points". Finite resistivity allows lines 
to reconnect at the nodal X-points and then contract toward 
the axes of the magnetic islands passing through the O-points. 
For a given wavenumber the instability can thus be visualized 
as an alternating set of perturbations in the current density, 6j, 
that tend to incre ase/decrease the local background curren t 
density (see, e.g., ISturrockl[T994t iBoyd & Sandersonll2003h . 
Although the hallmark of the tearing mode is its current den- 
sity pattern, the associated velocity field displays a distinct 
vorticity field pattern (see. Figure 5.7 in iBovd & Sanderson! 
12003.) . 

5.2. Vorticity and Current Density of Parasitic Modes 

In order to understand the physical nature of the secondary 
modes that are likely to play a role in the evolution of the 
primary MRI modes, it is of special interest to understand the 
fluid dynamics along the directions associated with the fastest 
growing parasites, 0,nax- As discussed in Section[3] for fixed 
values of the dissipation coefficients, the growth rates of the 
secondary instabilities peak around specific directions which 
are almost aligned with either the velocity or magnetic fields 



of the primary MRI mode, i.e., 9„ 



for A„ ^ 1 and 



9-Q for ^ 1, see Figure |3] Thus, motivated by 
the characteristics of the Kelvin-Helmholtz and tearing modes 
described above, it is of particular interest to understand the 
structure of the vorticity and current densities associated with 
the parasitic modes in the directions that are perpendicular to 
the planes defined by k\^^ and z, when 9 = (?max- 

In order to best illustrate the underlying physics of the 
fastest parasitic modes, it is useful to define a new coordinate 
system z) that is rotated with respect to the coordinate 
system defined by (z, y, z) with a rotation angle 6'„iax about 
the z direction. We define the coordinates h and p in such a 
way that they increase along the direction of the versors fch 
and fcp, respectively, see Figure|2] with 



feh X fep — z , 

Ajp X Z — ■ 



(75) 
(76) 



In this coordinate system, the velocity and magnetic fields of 
the parasitic modes have components 



5v ^ 5vhk\i + (5upfcp + 5vzZ , 
5B = SBhkh + SBpkp + SB^z . 



(77) 
(78) 



Since in these coordinates k^-x = k\^h, the secondary pertur- 
bations 5vz and SB^ in equations (ISTT l and ( |32] | become 



-i(n-\-kz)z —ikiih 



SB,{h,z)= J2 (3ne''^"+''''>''e"'"^^^ 



(79) 



(80) 



n— — oo 



^ In order to simplify notation, in this section we refer to tlie versor asso- 
ciated with the direction of fastest growth, i.e. fc^ niaxi simply as fc^ 



Because we are considering an incompressible fluid, we can 
calculate the projection of the velocity and magnetic field of a 
given parasitic mode onto the plane defined by (feh, z), i.e.. 



= {Sv-kh)kh + Sv^z , 
5Bu = {6B-ki,)ki, + SB,z, 



(81) 
(82) 



just in terms of the Fourier coefficients {an} and {/3„} ob- 
tained by solving the eigenvalue problem defined by Equa- 
tions ( |37] | and dJST l. Using that Sv-k}^ — dzSvz/iiki^) and 
SB-k\i — dzSBz/ {ikh) we obtain 



(83) 



OO 



n— — OO 



n^ — oo 



Note that the components of the three-dimensional velocity 
and magnetic fields that are orthogonal to these planes are 
independent of the coordinate p. Thus, using the fact that the 
divergence is invariant under rotations, we conclude that the 
divergence of the two-dimensional vector fields that lie on the 
plane (fch, z) should vanish, i.e., V-(5t>|| = V'(5J5|| = 0. 

The components of the vorticity and current density perpen- 
dicular to the plane defined by (fch, z), i.e., Su)± = Soj^kp 
and Sj = Sj±kp, are given by 

6luj_ = 5u3-kp = (V X 6v)-kp = ( V x 5v\\)-kp , (85) 

Sj±=Sj-kp - (V X SBykp = (V X ,5B||).fcp, (86) 

where the action of the curl operators in the rightmost equali- 
ties is given by 

V X 5v\\ = {dz5vh - dh5vz)kp , (87) 

V xSj^^^{dz6Bh~dh6Bz)kp. (88) 

Therefore, using equations ( [83] l and ( [84l l. we can obtain ex- 
pressions for Su)± and Sj ^ in terms of the Fourier coefficients 

{an} and {/?„}, 



ik 



ik 



(89) 



— [kl + {n + kzf](3ne- 



i(n-\-kz)z ^—ik^^h 



(90) 



5.3. Lagrangian Displacement induced by Parasitic Modes 

The vertical Lagrangian displacement provides useful 
complementary information to help us identify the physical 
nature of the eigenmodes. Let us then briefly describe how to 
calculate in terms of the Fourier coefficients {a„} that re- 
sult from solving the eigenvalue problem. The rate of change 
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k„h 

Fig. 8. — The arrows in the upper and lower panels correspond to the projections of the velocity, (5i)|| {h, z), and magnetic field, (5B|| {h, z), of the parasitic 
modes onto the plane defined by the z-axis and the direction Smax along which they exhibit their fastest growth. The color contours correspond to the associated 
vorticity, 5u}j_, and current density, Sjj^, projected in the direction perpendicular to 9max, i.e., along fcp, see Figs.[T]and|2] The red and blue colors con'espond 
to the maximum positive and minimum negative values (with the corresponding vectors pointing into and out of the page) associated with the vorticity and cuirent 
density of each mode. Here, A^/ = 10'^ and Pm = {10^^, 10^'^, 10^^, 10^^}, from left to right. The left-most panels, which correspond to A„ = 0.1, show 
the cuiTent density and vorticity patterns (lower and upper panels) characteristic of the tearing mode instability, (see. Fig. 5.7 in Boyd & Sanderson 2003). Note 
the convergence of the flow velocity toward the null surfaces of the MRI magnetic field, AB oc cos{Kz), at z = n-n 12. These type of motions are ultimately 
responsible for the reconnection of the magnetic field of the primary MRI mode and would eventually limit its growth, see also Fig.|9] The next set of three 
(upper and lower) panels correspond to Arj = {1, 10, 10^}, from left to right. The upper panels show clear signatures of the Kelvin-Helmholtz instability in 
the vorticity contours, see e.g., (see Fig. 1.7.2 in Batcheloi 2000). The (less familiar) lower panels show the cun'ent density perturbations associated with these 
Kelvin-Helmholtz modes. We note that the amplitudes of the fluctuations in the current density are much smaller than the ones in the vorticity. In this cases, the 
lower panels do not show fluctuations in the current density contours in the locations of the null magnetic surfaces. Thus, this (weak) current density pattern does 
not seem to lead to the reconnection of the MRI field. This is confirmed by the behavior of the vertical Lagrangian displacement, see Fig. [9] 



in the Lagrangian displacement £ with respect to a point mov- 
ing with velocity v is given by 

§ = ^ + (t;.V)| = ,5t; + ^.V7;. (91) 
at at 

In the periodic background provided by the primary MRI 
mode, the Lagrangian displacement is of the form ^(a;, t) = 
^o(z) exp [st — zfe-cc] with ^Q(cc + 27r/i4r2;) = ^o(^)- There- 
fore, considering the MRI velocity field, — sin(i^z), 
as the background velocity, we can obtain in terms of the 
vertical Eulerian velocity as 

Sv,(h,z) 

Uh,z) = 4r^- (92) 

5.4. Parasitic Mode Identification 

In Section |4] we showed that the eigenvalues correspond- 
ing to the growth rates of the parasitic modes reached well 
defined asymptotic regimes for ^ 1 and ^ 1 in 
the limit A,y ^ 1. This must also be true for the eigen- 
modes. It is then only necessary to explore in detail the be- 
havior of the modes for moderate values of A,,. We thus fo- 
cus our attention in the region of parameter space spanned by 
A^ = {10^^, 1, 10, 10^}. Consider as an example^ the set of 

* Because the scalings presented in Section |4] if we were to chose 
K = 10^ or A^ = 10'' we would obtain the same results for the range 



of magnetic Prandtl numbers Pm = {10 ,10 

{10~5,10~**,10~^,10"2| 



,10- 



, 1} or Pm 



parameters A^ = 10^ and Pm = {10"'', 10"^^ IQ-i}. 
The structure of the corresponding fastest parasitic modes are 
shown in Figure [8] from left to tight. The projections of 
the velocity and magnetic fields onto the planes defined by 
(fch, z), i.e., Sv\^ {h, z) and 5B\\ {h, z), are shown with white 
arrows in the upper and lower set of panels, respectively. The 
color contours show the projection of the vorticity and current 
density along the direction fcp (perpendicular to the page), i.e., 
Su)± and Sj j^. 

The versor characterizing the direction of fastest growth, 
fch, for the Elsasser number A,, = 10^^ points in the direction 
0,„ax — Sb, see Fig. (O and Table[T] This mode feeds off the 
current density of the primary MRI mode. The correspond- 
ing mode structure resembles closely the perturbations in the 
current density and induced vorticity patterns expected in the 
analysis of the stability of a set of equidistant current sheets 
distributed along the z direction and alternating sense accord- 
ing to ±fcp, (see. Figure 5.7 in Boyd & Sanderson 2003). The 
current density of the secondary modes presents maxima and 
minima along the planes z = ±rm/2 where the magnetic 
field of the primary mode, AS = Bq cos{Kz), reverses 
sign. Thus, the fluctuations induced by these fastest resis- 
tive secondary modes tend to promote reconnection of the 
MRI field^. This can be better appreciated by analyzing the 

^ There is another set of periodic maxima and minima in the current density 
fluctuations that lie on planes z = 0, itnTr; these are the locations where the 
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riff 







Fig. 9. — Vertical Lagrangian displacement ^z{h,z) coiresponding to the fastest growing parasitic instabilities projected onto the plane defined by the 
z-axis and the direction associated with fastest growth, Omax- The panels correspond to = 10^ and Pm = {lO"*, IQ-^, 10"^, 10"^}, i.e.. A,, = 
{10~^, 1, 10, 10'^}, from left to right. The convergence of the vertical displacement field toward the null surfaces of the MRI magnetic field, AB oc cos{Kz), 
located at z = itn7r/2, is evident in the leftmost panel. This feature is characteristic of tearing modes; the associated motions are ultimately responsible for 
reconnecting the MRI magnetic field and limiting its growth. The next set of three panels show the vertical displacements associated with the fastest Kelvin- 
Helmholtz modes, which mainly tend to bend the horizontal MRI magnetic field without necessarily promoting its reconnection (see Fig.[8]) 

sign. The corresponding vertical Lagrangian displacement as- 
sociated with these modes for = {1, 10, 10^} are shown in 
the three rightmost panels in Figure |9] respectively. Thus, the 
fluctuations induced by the fastest growing Kelvin-Helmholtz 
secondary instabilities tend to bend the horizontal MRI mag- 
netic field without directly promoting their reconnection. 

The upper and lower panels of Figure [TO] show the projec- 
tions of the total (primary plus secondary) velocity and mag- 
netic fields onto the planes defined by (fch, z), i.e., AV|| + 
Sv\\ [h, z) and ASy + 5B\\ (h, z), where 



vertical Lagrangian displacement defined in Equation 
which is shown in the leftmost panel of Figure |9] The ob- 
served mode structure is qualitatively insensitive to the value 
of the Elsasser number as long as A.,, < 1 and Ki, 3> 1, the 
growth rates and lengthscales associated with each value of 
A^ change, of course, as discussed in Section |4] We thus 
conclude that the fastest parasitic modes correspond to tear- 
ing modes for A^ < 1. These parasitic modes are enabled 
by non-zero resist ivity and are thu s absen t in the ideal MHD 
regime studied bv lGoodman & Xul (I1994I) . 

The versors fch characterizing the direction of fastest 
growth for the Elsasser numbers A,, = {1, 10, 10^} point in 
the direction 6'max — ^'v for the three rightmost panels in 
Fig. (O, see also Table [T] These modes feed off the shear in 
the velocity field of the corresponding primary MRI modes. 
The velocity and vorticity fields show a periodic structure 
similar to what is expected from the stability analysis of a pe- 
riodic set of equidistant vortex sheets distributed along the z 
direction and alternating sense according to ifcp. The mode 
structure in the velocity and vorticity fields, as well as the 
growth rates and lengthscales, associated with each value of 
A^ are quantitatively insensitive to the value of the Elsasser 
number as long as A.,, > 1 and A,y ^ 1, as discussed in Sec- 
tion g] 

The white arrows in the lower panels show the mag- 
netic field projected onto the plane defined by (fch, z). The 
color contours correspond to the associated current densities 
along the directions fcp. The structures seen in these pan- 
els show some evolution as a function of A^ > 1. How- 
ever, these modes are not intrinsically modified from pure 
Kelvin-Helmholtz modes since the amplitude of the fluctu- 
ations in the magnetic and current density fields is much 
smaller than the fluctuations in the velocity and vorticity 
fields. We thus conclude that the fastest parasitic modes cor- 
respond to Kelvin-Helmholtz modes for A^ > 1. In the limit 
A^ ^ 1, these corres pond of course to the K elvin-Helmholtz 
modes alluded to in iGoodman & Xul (119941) . Note that the 
current density of the secondary modes vanishes along, and 
also in the vicinity of, the planes z — ±n7r/2 where the mag- 
netic field of the primary mode, AS = Bq cos{Kz), changes 

currents associated with the MRI magnetic field vanish. Thus, the presence of 
these current density perturbations seem to be needed to satisfy the periodic 
constraints (on the scale of the unstable mode) rather than due to the unstable 
configuration presented by the MRI cun'ents themselves. 



Ay II = Vo sin(Xz) cos(0,nax - 0v) , (93) 
AB|| = Bo cos(i^z) cos(0„,ax - 0b) . (94) 

The four panels correspond to A^ = {10^^, 1, 10, 10^}, from 
left to right. The color contours show the projection of the 
total vorticity and current density along the direction fcp (per- 
pendicular to the page), i.e., 6u)±fl + 5u3± and 5j q + 6j ^. 
The contributions to the vorticity and current density associ- 
ated with the background velocity and magnetic fields of the 
primary MRI mode are given by 

Sujo.x = V X AF|| = Vo cos{Kz) co8{e ~ 6^) fcp , (95) 
<5jo,± = V X AB|| = -Bq sin(ifz) cos(6i - 6^) fcp (96) 

In all cases, the amplitude of the primary mode Vo and Bq 
is such that the fastest secondary modes grow as fast as the 
primary MRI mode, i.e., s^^^iy, 77, A'^ax) = rf). 

The addition of the MRI fields facilitates the identification 
of the velocity and magnetic fields that result from the in- 
fluence of the secondary instabilities with the more familiar 
structures that are expected from Kelvin-Helmholtz and tear- 
ing mode instabilities in periodic backgrounds. In particular, 
it highlights the presence of "O" and "X" points in the case 
of the tearing mode and the wave-like structure of the veloc- 
ity field in the cases associated with Kelvin-Helmholtz modes. 
Note that the contributions to the total velocity and vorticity 
field from the MRI mode is relatively small in the case of 
= 10^^ (i.e., the tearing mode in the leftmost panel). This 
is because the ratio Vo/i?o is small and O^^^x ~ Oy ~ 7r/2. 
Analogously, the contributions to the total magnetic and cur- 
rent density field from the MRI mode is relatively small in 
the case of A^ — 10^ (i.e., the Kelvin-Helmholtz mode in the 
rightmost panel). In this case, 6'inax — ^'b — — 7r/2. 
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Fig. 10. — The arrows in the upper and lower panels correspond, respectively, to the projections of the total (primary plus secondary) velocity, 
AV|| + (h, z), and magnetic fields, ABy + (h, z), onto the plane defined by the z-axis and the direction 9max associated with the fastest par- 
asitic modes, see Eqs. j93) and i94\ . The amplitudes of the MRI velocity Vq and magnetic fields Bg correspond to the values at which the fastest sec- 
ondary modes reach a growth rate equal to the fastest growing MRI mode, for the con'esponding values of the dissipation coefficients. Here, = 10'^ 
and Pm = {10"'*, 10~^, 10~^, 10~^}, i.e.. A,, = {10^^, 1, 10, 10^}, from left to right. The color contours correspond to the associated total vorticity, 
5u}±fi + Su!±, and cuirent density, 5jj^ q -f Sj±, see Eqs. (95) and i96\ . projected in the direction perpendicular to 9niax, i.e., along the versor fcp, see Fig.|2] 
The red and blue colors correspond to the maximum positive and minimum negative values (with the coiresponding vectors pointing into and out of the page) 
of the total vorticity and total current density associated with each parasitic mode. The addition of the MRI magnetic field makes the presence of structures 
that resemble the "O" and "X" points, characteristics of the tearing mode instability, more evident in the leftmost, lower panel for which A^ = 1 (see, e.g. 
[Sovd & Sanderson 2003). Analogously, the inclusion of the MRI velocity field strengthens the identification of the modes coiresponding to A,j = {1, 10, 10^} 
with Kelvin-Helmholtz instabilities. 



6. SUMMARY AND DISCUSSION 

We have solved for the spectrum of parasitic instabilities 
that feed off exact MRI modes for a broad range of dissipa- 
tion coefficients. We parameterized the effects of viscosity 
and resistivity via the Elsasser number number A,, and its vis- 
cous counterpart A^. Our approach to this problem enabled 
us to explore dissipative regimes that are relevant to astro- 
physical and laboratory conditions that lie beyond the regime 
accessible to current numerical simulations. We explored in a 
systematic way the physical properties of parasitic modes for 
a grid of parameters defined by = {1, 10, ... , 10^} and 
magnetic Prandtl numbers Pm = {10^^, 10~^, . . . , 10^}. 
This allowed us to identify the existence of a critical Elsasser 
number of order unity and two distinct asymptotic regimes for 
quasi-ideal, ^ 1, and inviscid, resistive MHD, A^ ^ 1. 

The dynamics associated with arbitrary combinations of the 
various parameters involved (six after choosing the fastest 
MRI mode to set the growth rate and g eometry of the pri - 
mary), can be in general rather complex (iLatter et al.ll2009t) . 
We have thus chosen to focus our attention on the fastest 
growing parasitic modes. Axisymmetric modes, i.e., 9 = 0, 
present growth rates that can be rather small compared to the 
fastest parasites. However, they might deserve some atten- 
tion as they seem to play an imp ortant role in the saturation of 
2D (axisymmetric) simulations (ISanoll2007l; iMasada & Sand 
120081: lObergaulinger et al.ll2009l) . 

We showed that important properties of the fastest sec- 



ondary instabilities found in iPessah & GoodmanI (l2009h for 
^ ^ K < 100 and 0.1 < Pm < 10 are generic. The 
fastest parasitic modes are non-axisymmetric, have purely 
real growth rates, have the same vertical periodicity as the pri- 
mary mode upon which the feed, and horizontal wavelengths 
that are roughly twice as long. Their wavevectors fch are 
nearly aligned with either the velocity or the magnetic field 
of the primary mode. The first type dominate for A,, > 1 
and correspond to Kelvin-Helmholtz modes that feed off the 
vertical shear induced by the MRI and reach their maximum 
growth rates along the direction of the MRI velocity field. The 
second type dominate for A,, < 1 and are related to tearing 
modes that feed off the MRI currents and grow the fastest 
along the direction of the MRI magnetic field. We provided a 
set of scaling laws that describe the growth rates of the fastest 
growing parasitic Kelvin-Helmholtz and tearing instabilities 
in two important asymptotic regimes 



no 
7^7 



:0.20 



Bo 



for A,, > 1 



A,, > 



1, 



Bo . 
^0.15-g^A, 



for Ai. > 1 , A„ < 1 



(97) 



(98) 



Following iPessah & GoodmanI (l2009t) . we consider the 
"saturation amplitude" B^^{v,ri) as the amplitude that the 
magnetic field produced by the MRI must have grown to in 
order for the instantaneous growth rate of the fastest parasitic 
mode, Snia.y^{v, rj), to match that of the fastest primary mode, 
^maxi'^,'!])- Because both the growth rate of the parasitic 
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modes and the MRI decay proportionally to for A^ <C 1, 
the amplitude to which the MRI magnetic field needs to grow 
to in order for these growth rates to become comparable is 
independent of the Elsasser number As resistivity increases 
the magnetic field that can be generated by the MRI before 
the tearing modes become dynamically important reaches an 
asymptotic value. However, due to the alignment of the MRI 
magnetic field with the azimuthal direction at low Elsasser 
number, the situation is quite different for the stresses. The 
stresses produced by the MRI primary modes, at the time 
when the fastest tearing modes have growth rates similar to 
their own, decay linearly with decreasing Elsasser number. 



i.e., Tl 



(X A^ for A,, 



^ 1. This behavior seems consis- 



tent with numerical simulations of resistive MHD shearing 
boxes, where the stresses at sat uration decrea se linearly with 
Elsasser number (see Table 1 in lFleming et al.i.2000: note that 
the associated Elsasser numbers in the more resistive runs are 
somewhat larger than unity). 

In the region of parameter space where our analysis over- 
laps with the regime accessible to numerical simulations, 
the values of asat are within factors of a few of the val- 
ues obtained in the turb ulent regime. As pointed out in 
iPessah & Goodm an (2009). the dependence on the stresses at 
saturation on the magnetic Prandtl number is less pronounced 
and the predicted value of asat is smaller, by a factor of 6 
at Pm = 1, tha n that found in the non-linear simulations of 
iLesur & Longare tti (2007). Part of these differences can be 
accounted for with a mor e sensible operational definition of 
saturation. iPessah & Goodman (.2009) provide a way to es- 
timate the amplitude of the fields at saturation when primary 
and secondary modes reach comparable amplitudes based on 
the values obtained when they reach equal growth rates. They 
conclude that the overshoot factor, which depends logarithmi- 
cally on the initial amplitude of the parasite, is likely between 
~ 3 and ~ 10. These arguments should also be applicable at 
low Elsasser numbers since they mainly rely on the fact that 
the growth rate of the secondary modes is linear in the am- 
plitude of the primary MRI mode. It is then conceivable that 
the asymptotic analysis carried out in this paper could pro- 
vide, within factors of a few, a reasonable description of the 
saturation level in a wide region of parameter space. 

We explored in detail the physics of the eigenmodes with 
fastest growth rates in the regime where variations in the El- 
sasser number lead to non-trivial modifications to the struc- 
ture of the parasitic modes. To this end we emphasized the 
importance of understanding the structure of the vorticity and 
current density patterns associated with the secondary insta- 
bilities and used these as a mean to confirm their association 
with Kelvin-Helmholtz and tearing mode instabilities. Our 
analysis suggests a strategy to recognize the presence of par- 
asitic modes and attempt to identify their nature: (i) Evolve 
the simulation until the breakdown of the initial exponential 
growth or subsequent peaks in stress or magnetic energy, (ii) 
Dissect the simulation domain in planes perpendicular to the 
mid-plane that contain the z-axis. (iii) Project the velocity 
and magnetic field in these planes and take the curls in order 



to obtain the corresponding vorticity and current density, (iv) 
Determine whether these resemble what is expected from the 
Kelvin-Helmholtz or tearing mode instabilities discussed in 
Section|5] 

Understanding from first principles the saturation of the 
MHD turbulence, and the associated angular momentum 
transport, driven by the MRI under realistic physical condi- 
tions (including, ambipolar diffusion, compressibility, strati- 
fication, radiation, etc.) is a challenging endeavor The re- 
sults presented throughout this study provide guidance to- 
ward understanding some of the key physical processes that 
are likely to play a role in halting the exponential growth of 
the M RI and, mo re speculative ly, in the subsequent turbulent 



state (S ano & Inutsuka 2001,; iSanol 



2007tlBodoetal.| l2008t 
IPessah & Goodman 2009, see also iLatter et al .1120091) . For 
instance, there is numerical evidence suggesting that ohmic 
heating due to the reconnection of MRI field lines is an im- 
portant source o f energy in shearing box simulations of resis- 
tive MHD (Fleming et aT 2000; Sano & Inutsuka 2001). In 
allowing for non-zero resistivity and calculating the currents 
associated with the secondary instabilities, we have gone one 
step forward in establishing the chain of processes that enable 
the conversion of gravitational energy into thermal energy in 
differentially rotating, magnetized, non-ideal plasmas. 

As a final comment, we mention two limitations that re- 
sult from our working assumptions that are worth emphasiz- 
ing due to their potential si gnificance (see al so the related 
discussion in Section 2.4 in iLatter et ar]|2009l) : (i) We have 
neglected the effects of shear on the dynamics of the sec- 
ondary modes. Non-axisymmetric para sitic modes will shear 
linearly in time (iGoodman & Xuiri994l) : therefore, assuming 
fixed horizontal versors to characterize a parasitic mode is 
an approximation, (ii) We have not considered the explicit 
coupling between the evolution of the MRI-modes and the 
secondary instabilities. As the secondary modes grow they 
drain energy from the primary modes; therefore the secondary 
growth rates, that rely on the amplitude of the MRI modes, 
could be affected. Both of these effects are likely to diminish 
the ability of the parasitic modes discussed here to slow the 
growth of the MRI, or the subsequent "channel flows". Thus 
the extrapolation of the results presented here to the non-linear 
regime should be complemented with the pertinent quota of 
skepticism. Having said this, at present, the properties of par- 
asitic modes provide valuable analytical guidance and a basic 
framework to design and interpret tailored numerical experi- 
ments in order to shed light into the non-linear saturation of 
the MRI. 
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